a decay half-lives of new superheavy nuclei within a generalized liquid drop model 
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The a decay half-lives of the recently produced isotopes of the 112, 114, 116 and 118 nuclei and 
decay products have been calculated in the quasi-molecular shape path using the experimental Q a 
value and a Generalized Liquid Drop Model including the proximity effects between nucleons in the 
neck or the gap between the nascent fragments. Reasonable estimates are obtained for the observed 
a decay half-lives. The results are compared with calculations using the Density-Dependent M3Y 
effective interaction and the Viola-Seaborg-Sobiczewski formulae. Generalized Liquid Drop Model 
predictions are provided for the a decay half-lives of other superheavy nuclei using the Finite Range 
Droplet Model Q a and compared with the values derived from the VSS formulae. 
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The synthesis of superheavy elements has advanced us- 
ingboth cold [l| and warm fusion reactions. Recently 
llBlil, isotopes of the elements 112, 114, 116 and 118 
have been produced in fusion-evaporation reactions at 
low excitation energies by irradiations of the 233 > 238 u, 
242 Pu, 248 Cm and 249 Cf targets with 48 Ca beams. The 
main decay mode is the a emission and the a decay en- 
ergies and half-lives of fourteen new a decaying nuclei 
have been measured. Some questions have been raised 
Q about these superheavy element findings. In similar 
sophisticated experiments at other places 00 the a cas- 
cades were not observed. 

The pure Coulomb barrier sharply peaked at the touch- 
ing point does not allow to determine correctly the par- 
tial a decay half-lives. It is probable that the a decay 
takes place in the quasi-molecular shape path where the 
nucleon-nucleon forces act strongly during the formation 
of the neck between the nascent fragments and after the 
separation and a proximity energy term must be added in 
the usual development of the liquid-drop model Q . The 
generalized liquid drop model (GLDM) which includes 
such a proximity e nerg y term has allowed to describe the 
fusion j| , fission [lOj . light nucleus an d o. emission 
processes. 

The purpose of this work is to determine the partial a 
decay half-lives of these superheavy elements within this 
GLDM from the experimental Q a values using the WKB 
approximation and to compare with the experimental 
data and the calculations with the Density-Dependent 
M3Y (DDM3Y) effective interaction Q and the Viola- 
Seaborg formulae with Sobiczewski constants (VSS) . 
Finally predictions within the GLDM and VSS formulae 
are given for the partial a decay half-lives of the still non 
observed superheavy nuclei ranging from Sg to Z = 120. 

For a deformed nucleus, the macroscopic GLDM en- 



ergy is defined as 

E = Ey + Es + Ec + -ERot + -Bprox- (1) 

When the nuclei are separated: 

E v = -15.494 [(1 - l.MDAx + (1 - 1.8/|)A 2 ] MeV, 

(2) 

MeV, 
(3) 



E s = 17.9439 



:i - 2.6/ 2 )A 2/3 + (1 - 2.6/f K 



E c = 0.6e 2 Z(/Ri + 0.6e 2 Z%/R 2 + e l Z 1 Z 2 /r, (4) 

where A,-, Zi, Ri and Ii are the mass number, charge 
number, radii and relative neutron excesses of the two 
nuclei, r is the distance between the mass centres. The 
radii Ri are given by : 



Ri = (1.28AV 3 - 0.76 + 0.8A,- v J ) fm 



-l/3\ 



(5) 



This formula allows to follow the experimentally observed 

1/3 

increase of the ratio rj=Ri/A i with the mass; for ex- 
ample, r = 1.13 fm for 48 Ca and r = 1.18 fm for 248 Cm. 

For one-body shapes, the surface and Coulomb ener- 
gies are defined as: 

E s = 17.9439(1 - 2.6/ 2 )A 2/3 (S'/47r J R 2 ) MeV, (6) 



E c = 0.6e 2 (Z 2 /i? )x0.5 / (V{9)/V )(R(6)/R ) 3 sin 9d0. 
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(7) 

S is the surface of the one-body deformed nucleus. V(8) 
is the electrostatic potential at the surface and Vq the 
surface potential of the sphere. 

The rotational energy is determined within the rigid- 
body ansatz: £R, t = h 21 ±^ • surface energy results 
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from the effects of the surface tension forces in a half 
space. When there are nucleons in regard in a neck or 
a gap between separated fragments an additional term 
called proximity energy must be added to take into ac- 
count the effects of the nuclear forces between the close 
surfaces. This term is essential to describe smoothly the 
one-body to two-body transition and to obtain reason- 
able fusion barrier heights. It moves the barrier top 
to an external position and strongly decreases the pure 
Coulomb barrier. 

^ProxM = 2 7 / $ [D(r, h)/b] 2nhdh, (8) 
A min 

where h is the distance varying from the neck radius or 
zero to the height of the neck border. D is the distance 
between the surfaces in regard and b — 0.99 fm the sur- 
face width. $ is the proximity function of Feldmeier |16| . 
The surface parameter 7 is the geometric mean between 
the surface parameters of the two nuclei or fragments. 
The combination of the GLDM and of a quasi-molecular 
shape sequence has allowed to reproduce the fusion bar- 
rier heights and radii, the fission and the a and cluster 
radioactivity data. 

For the a emission this very accurate formula simulates 
the proximity energy |l7^ : 

£ prox (r) = (4^ 7 )e- 1 - 38 ( r - fl -- fld )[0.6584A 2 / 3 (9) 

_(2_l!H + o.4692A 1 /3)r 
A 1 ' 3 

-0.02548A 1/3 r 2 + 0.01762r 3 ]. 

To obtain the a decay barrier from the contact point 
between the nascent a particle and daughter nucleus it 
is sufficient to add this proximity energy to the Coulomb 
repulsion. 

The half-life of a parent nucleus decaying via a emis- 
sion is calculated using the WKB barrier penetration 
probability. In a unified fission model, the decay con- 
stant of the a emitter is simply defined as A = vqP. The 
assault frequency vq has been taken as vq = 10 20 s _1 . 
The barrier penetrability P is calculated within the ac- 
tion integral 

P = exp[-- ^2B(r)(E{r) - E(sphere))). (10) 

" JRin 

The deformation energy (relative to sphere) is small until 
the rupture point between the fragments [T3| and the two 
following approximations may been used:i?i n = Rd + Ra 
and B(r) = fi where \i is the reduced mass. R ou t is sim- 
ply e 2 ZdZ Q /Q Q . The partial half-life is related to the 
decay constant A by Ti/ 2 = ^ ■ The a decay half- 
lives of the recently produced superheavy nuclei calcu- 
lated with the GLDM using the experimental Q a value 
and without considering the rotational contribution are 
presented in Table 1. The results agree reasonably with 
the experimental data indicating that a GLDM taking 
account the proximity effects, the mass asymmetry, and 



an accurate nuclear radius is sufficient to reproduce the a 
decay potential barriers when the experimental Q a value 
is known. The results obtained with the DDM3Y inter- 
action agree with the experimental data as the GLDM 
predictions and largely better than the VSS calcula- 
tions. This shows that a double folding potential ob- 
tained using M3Y 0] effective interaction supplemented 
by a zero-range potential for the single-nucleon exchange 
is very appropriate because its microscopic nature in- 
cludes many nuclear features, in particular a potential 
energy surface is inherently embedded in this descrip- 
tion. This double agreement shows that the experimental 
data themselves seem to be consistent. For most nuclei 
the predictions of the VSS model largely overestimate 
the half lives. The blocking effect is probably treated too 
roughly. 

The half live of 294 118 is slightly underestimated in the 
three theoretical calculations possibly due to the neutron 
submagic number N = 176. In Ref. yjj, it is also pointed 
out that for oblate deformed chain of Z= 112, the shell 
closure appears at N= 176. 

Most of the theoretical half lives using GLDM are 
slightly smaller than the experimental data. A reason 
is perhaps that the rotation of the nuclei is neglected in 
the present calculations. The term h 2 l(l + 1)/(2I^) in 
Eq.(??) represents an additional centrifugal contribution 
to the barrier which reduces the tunnelling probability 
and increases the half lives. A second reason is that the 
shell effects and pairing correlation are not explicitly in- 
cluded in the alpha decay barrier, in spite of their global 
inclusion in the decay energy Q. 

The experimental a decay half-lives are between the 
close theoretical values given by the GLDM and the ones 
derived from the VSS formulae. Thus predictions of the 
a decay half lives with the GLDM and VSS formulae are 
possible. In the next calculations the experiential Q a 
values are taken from the FRDM [2(j which reproduces 
all known experimental data of ground state properties of 
a large number of nuclei and gives good predictions for 
nuclei far from the (3 stability line and the superheavy 
nucleus region. In Ref. T1/2 obtained using the Q a 
values given by the Thomas-Fermi model can be found. 

a decay half-lives for Z — 106 to Z = 120 isotopes are 
shown in FigQ] the open dots indicating the results of 
GLDM and the black triangles the ones derived from the 
VSS formulae. 

The half-life values vary from years to microseconds. 
A narrow window exists between the two predictions for 
each isotope and the unknown a decay half-lives of SHN 
may lie in this window assuming that the FRDM Q a is 
correct. 

The FRDM Q a value explicitly displays a minimum at 
the submagic number N= 162. It is progressively eroded 
by the neutron deficiency and the subshell disappears 
completely from Z= 115. This induces a small first hump 
in the predicted logiopl^s)] curves till Z= 114. Before 
this first submagic number the half-lives increases rapidly 
with A from Db to Z— 114 isotopes. For the results of 
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TABLE I: Comparison between experimental a decay half-lives 4j and results obtained with the GLDM, the DDM3Y effective 
interaction [14} and the VSS formulae. 



Parent 


Nuclei 


Expt. 


Expt. 


DDM3Y 


GLDM 


VSS 


Z 


A 


Q(MeV) 


Tl/2 


Tl/2 


7l/2 


T1/2 


118 


294 


11.81 ± 0.06 


1.81^3 ms 


°- 66 -ai8 ms 


0.15+.^ ms 


rrmrra 

°- 64 -ai8 ms 


116 


293 


10.67 ±0.06 


53l^g ms 


^ r , -X- on 

206Ig? ms 


„ _ i in 99 

22.8lt7° 6 ms 


1258±§|4 ms 


116 


292 


10.80 ±0.07 


18le 6 ms 


_ „ i o n 

39li° ms 


IO.45I3 45 ms 


49lie ms 


116 


291 


10.89 ±0.07 


Q.Stfs ms 


_ „ , ion 9 

60.4l|°;i ms 


6.35j£]£ ms 


336.4+"^ ms 


116 


290 


11.00 ± 0.08 


15+f ms 


, „ . _|_7 7 

13.4+52 ms 


„ , 1 nil 

3.47±i;l| ms 


15.2+5;g ms 


114 


289 


9.96 ± 0.06 


^•'-0.7 S 


3.8_j_;| s 


„ p-rj+0.25 

U.O^ q ^7 s 


26.7lg 3 7 1 s 


114 


288 


10.09 ±0.07 


0.8+g-g s 


0.67 + E! : ^ s 


0.22 q Qg s 


0.98 + giS s 


114 


287 


10.16 ±0.06 


0.511°;!* s 


1 n+ - 52 s 

1 - lo -0.40 b 


0.161S& s 


7.24±*£ s 


114 


286 


10.35 ±0.06 


o.i6±g:8l s 


0.14i° ;°« s 


0.05±°;g| s 


0.19+11 s 


112 


285 


9.29 ± 0.06 


34t* 7 s 


75+« s 


13.2211'^ s 


592 2Q7 ^ 


112 


283 


9.67 ±0.06 


4.0lJ; 3 7 s 


5.9+11 s 


0.95+SH s 


41.3+™ s 


110 


279 


9.84 ± 0.06 


18+ - 05 s 


0.40+.^ s 


o.oal°£ s 


2 Q2+ 1 - 4 s 


108 


275 


9.44 ± 0.07 


u - io -0.06 H 


1 nq+ 73 q 
i.uy_ 040 s 


o.27±S;!S s 


8.98t|;l s 


106 


271 


8.65 ± 0.08 


2A±f 3 Q min 


1.0±g;| min 


0.33±g;?| min 


8.6+.H min 
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FIG. 1: Comparison between calculated a-decay half-lives of Sg, Hs, Ds, 112, 114, 116, 118 and 120 isotopes using the GLDM 
and the VSS formulae. 



VSS formulae, it seems that the turning point is delayed 
a little and appears at neutron number N= 163, for the 
blocking effect has been magnified and the half lives are 
overestimated. N= 184 is always a closure shell. Subclo- 
sure or closure shells exist also around N= 176 for Z= 105 
to Z= 112. The most stable nuclei should stand about 
N= 184 when the proton number is higher than Z= 115, 
but stand around N= 176 for lower Z values. 

As a conclusion the half-lives for a-radioactivity have 
been analyzed in the quasimolecular shape path within a 
Generalized Liquid Drop Model including the proximity 
effects between nucleons and the mass and charge asym- 
metry. The results are in reasonable agreement with the 



published experimental data for the alpha decay half- 
lives of isotopes of charge 112, 114, 116 and 118 and close 
to the ones derived from the DDM3Y effective interac- 
tion. The experimental a decay half lives stand between 
the GLDM calculations and VSS formulae results and 
the a decay half-lives of still non-observed superheavy 
nuclei have been predicted within the GLDM and VSS 
approaches and Q Q derived from the FRDM. 
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